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Abstract 

According to a rederivation - due to Collins and Qiu - the DGLAP equation can 
be reinterpreted (in leading order) in a probabilistic way. This form of the equation 
has been used indirectly to prove the bound \Af(x, Q)\ < f(x, Q) between polarized 
and unpolarized distributions, or positivity of the helicity distributions, for any Q. 
We reanalize this issue by performing a detailed numerical study of the positivity 
bounds of the helicity distributions. To obtain the numerical solution we implement 
an x-space based algorithm for polarized and unpolarized distributions to next-to- 
leading order in a s , which we illustrate. We also elaborate on some of the formal 
properties of the Collins-Qiu form and comment on the underlying regularization, 
introduce a Kramers-Moyal expansion of the equation and briefly analize its Fokker- 
Planck approximation. These follow quite naturally once the master version is given. 
We illustrate this expansion both for the valence quark distribution qy and for the 
transverse spin distribution h\. 



QCD, the theory of the strong interactions, has reached a stage in which precision mea- 
surements of its dynamics have become possible. This scrutiny has allowed to obtain a 
better understanding of the fundamental structure of the nucleons, disentangling some 
important features of the underlying quark-gluon interaction with a very good accuracy. 

A lot of effort has been undertaken in the last few years to extend the same picture - 
at the same level of accuracy - also to polarized collisions, with a systematic perturbative 
analysis performed up to next-to-leading order (NLO) and, in part, in the unpolarized 
case, also to next-to-next to leading (NNLO) in a s , the QCD coupling constant. 

The aim of this theoretical and experimental effort, in particular at RHIC, the Rela- 
tivists Heavy Ion Collider at Brookhaven, is to describe the polarized spin distributions 
of the nucleon with accuracy. 

Therefore the study of possible theoretical constraints on the form of the initial con- 
ditions for these distributions and their evolution under the renormalization group (RG) 
turns out to be very useful. 

An interesting constraint relating longitudinally polarized, unpolarized and trans- 
versely polarized distributions is Soffer's inequality, which deserves a special attention, 
since has to be respected by the evolution to any order in a s . Some tests of the inequality 
have been performed in the near past, bringing support to it. However, other inequalities 
are supposed to hold as well. 

In this work we perform a NLO analysis of an inequality which relates longitudinally 
polarized distributions and unpolarized ones. The inequality can be summarized in the 
statement that helicity distributions (positive and negative) for quarks and gluons have 
to be positive. The inequality states that 

|A/(x,g 2 )|</(x,g 2 ) (i) 

or 

/±(x,Q 2 )>0 (2) 

where the ± refers to the the possible values of the helicities of quarks and gluons. The 
statement is supposed to hold, at least in leading order, for any Q. To analize the renor- 
malization group evolution of this relation, especially to next-to-leading order, requires 
some effort since this study involves a combined study of the (longitudinally) polarized 
and unpolarized evolutions. In this work we present a complete NLO study of the evo- 
lution equations starting directly from the helicity basis. Helicities are in fact the basic 
parton distributions from which other distributions can be built. 

Compared to other implementations, in our work we perform a NLO test of the pos- 
itivity of the helicity distributions using an ansatz due to Rossi [E] which reduces the 
evolution equations to an infinite set of recursion relations for some scale invariant coef- 
ficients. We have developed a complete implementation of this algorithm which will be 
made available and documented in related work of ours [15J. 

Various arguments to validate eq. (J2j) have been presented in the literature. From our 
perspective, an interesting one has been formulated by Teryaev and Collaborators who 
have tried to establish a link, to leading order, between evolution equations and their 
probabilistic interpretation in order to prove Soffer's inequality. Similar arguments hold 
also in the analysis of eq. ((2j). 
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oicier unpoiarizeo evolution |rj aim trie arguments 01 [Zj are inspired oy tne tact tnat 
the subtraction terms (the x — 1 contributions in the expressions of the kernels, where 
x is Bjorken's variable), being positive, once they are combined with the bulk (x < 1) 
contributions give a form of the evolution equations which are diagonal in parton type 
and resemble "kinetic" equations. Our arguments, on this issue, are just a refinement of 
this previous and influential analysis. 

In the recent literature there has been some attention to this feature of the DGLAP 
evolution, limited to the non-singlet sector, in connection with kinetic theory and the 
"dynamical renormalization group", in the words of ref. 0]. 

All the arguments, so far, go back to some important older work of Collins and Qiu who 
provided an interesting derivation of the (unpolarized) DGLAP equation using Mueller's 
formalism of cut diagrams. In their paper [T] the authors reinterpreted the DGLAP 
equation as a kinetic probabilistic equation of Boltzmann type. The authors gave no detail 
on some of the issues concerning the regularization of their diagrammatic expansion, on 
which we will elaborate since we need it for our accurate numerical analysis. In our work 
the Collins-Qiu form of the DGLAP equation is interpreted simply as a master equation 
rather than a Boltzmann equation, given the absence of a 2-to-2 scattering cross section 
in the probabilistic partonic interpretation. A master equation is governed by transition 
probabilities and various formal approximations find their way once this conceptual step is 
made. We illustrate, in the spirit of a stochastic approach to the DGLAP dynamics, how 
to extract standard differential equations of Kramers-Moyal type for the simplest non- 
singlet evolutions, those involving valence distributions and transverse spin distributions. 
Our analysis on this point is self-contained but purposely short, since a more detailed 
numerical and formal study of this developement is under way [T5] . 

We show that the DGLAP dynamics can be described, at least in a formal way, by a 
differential equation of arbitrarily high order. Truncations of this expansion to the first 
few orders provide the usual link with the Fokker-Planck approximation, the Langevin 
equation and its path integral version 1 . The picture one should have in mind, at least 
in this approximation, is that of a stochastic (brownian) dynamics of Bjorken's variable 
x in a fictititious time \og(Q), describing the evolution under the renormalization group 
(RG). In this interpretation the probability function is the parton distribution itself. 

2 Master Equations and Positivity 

Let's start considering a generic 1-D master equation for transition probabilities w(x\x') 
which we interpret as the probability of making a transition to a point x given a starting 
point x' for a given physical system. The picture we have in mind is that of a gas of 
particles making collisions in 1-D and entering the interval (x, x + dx) with a probability 
w(x\x') per single transition, or leaving it with a transition probability w(x'\x). In general 
one writes down a master equation 

d f 

—f(x 1 r) = J dx' (w(x\x')f(x', t) — w(x'\x)f(x, r)) dx'. (3) 

1 For an example of this interplay between differential and stochastic descriptions we refer the reader 
to 
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01 a many repiicas 01 waiKers 01 density j \x,T) jumping witn a pre-assigneu prouaomty, 
according to taste. 

The result of Collins a Qiu, who were after a derivation of the DGLAP equation 
that could include automatically also the "edge point" contributions (or x=l terms of the 
DGLAP kernels) is in pointing out the existence of a probabilistic picture of the DGLAP 
dynamics. These edge point terms had been always introduced in the past only by hand 
and serve to enforce the baryon number sum rule and the momentum sum rule as Q, the 
momentum scale, varies. 

The kinetic interpretation was used in [2] to provide an alternative proof of Soffer's 
inequality. We recall that this inequality 

\h 1 (x)\<q + (x) (4) 

famous by now, sets a bound on the transverse spin distribution hx(x) in terms of the 
components of the positive helicity component of the quarks, for a given flavour. The 
inequality has to be respected by the evolution. We recall that h±, also denoted by the 
symbol 

A T q(x,Q 2 )=q\x,Q 2 )-q l (x,Q 2 ), (5) 

has the property of being purely non-singlet and of appearing at leading twist. It is 
identifiable in transversely polarized hadron-hadron collisions and not in Deep Inelastic 
Scattering (DIS), where can appear only through an insertion of the electron mass in the 
unitarity graph of DIS. 

The connection between the Collins-Qiu form of the DGLAP equation and the master 
equation is established as follows. The DGLAP equation, in its original formulation is 
generically written as 

dq{x,Q 2 ) ( x dy 2 

where we are assuming a scalar form of the equation, such as in the non-singlet sector. 
The generalization to the singlet sector of the arguments given below is, of course, quite 
straightforward. To arrive at a probabilistic picture of the equation we start reinterpreting 
r = \og(Q 2 ) as a time variable, while the parton density q(x, r) lives in a one dimensional 
(Bjorken) x space. 

We recall that the kernels are defined as "plus" distributions. Conservation of baryon 
number, for instance, is enforced by the addition of edge-point contributions proportional 
to 5(1 - x). 

We start with the following form of the kernel 

Ptz) = P(z) - 5(1 - z) f 1 P(z) dz, (7) 

Jo 

where we have separated the edge point contributions from the rest of the kernel, here 
called P(z). This manipulation is understood in all the equations that follow. The 
equation is rewritten in the following form 

d ( \ f 1 r fi(x\q(y,r) f x dy ~ fy\ q(x,r) 
-q(x, r) = I dyP (- j — ^ -P (- j _ (8) 
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Figure 1: The constrained random walk of the parton densities 
Now, if we define 

/ I \ oi s ~ 6(y > x) 

w{x\y) = —P(xy) (9) 

2ti y 

(JBJ becomes a master equation for the probability function q(x,r) 
d f 

—q(x,r)= J dx (w(x\x)q(x, r) — w(x'\x)q(x, r)) dx . (10) 

There are some interesting features of this special master equation. Differently from other 
master equations, where transitions are allowed from a given x both toward y > x and 
y < x, in this case, transitions toward x take place only from values y > x and leave the 
momentum cell (x, x + dx) only toward smaller y values (see Fig.JT|). 

Clearly, this sets a direction of the kinetic evolution of the densities from large x values 
toward smaller-x values as r, the fictitious "time" variable, increases. 

Probably this is the simplest illustration of the fact that parton densities, at large 
final evolution scales, are dominated by their small-x behaviour. As the "randomly moving 
partons" reach the x ~ region of momentum space, they find no space where to go, while 
other partons tend to pile up toward the same region from above. This is the picture of a 
random walk biased to move downward (toward the small-x region) and is illustrated in 
Fig. ©. 



3 Probabilistic Kernels 

We briefly discuss some salient features of the structure of the kernels in this approach 
and comment on the type of regularization involved in order to define them appropriately. 
We recall that unpolarized and polarized kernels, in leading order, are given by 

PHI = Pff = C F ( JT ^ r -l-x + ls{l-x)) 



where 



and 



P<J) = 27)(z 2 + (l- 
p (0) = c y l + (!-*)' 



m 



x 



Cf = ^4, T f = T R n f = \n h (3 = ^N c - ±T f (12) 



Ap(°) - AP(0) 
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APW = 2T / (2x-l) 
APW = 0(2 -x) 

APW = 2iV c ( ir ^--2x + l) + 5(l-x)|, (13) 



while the LO transverse kernels are given by 



^ p «= c 4jrhr + - 2+ l S{l - x) )- (14) 



The unpolarized kernels should be compared with the Collins-Qiu form 

Pqq = Iqq-SQ-X) f dz^ qq 

J 

p gg = laa ~ ( n / J Q dz lqg + \ J Q dz lgg^j K 1 ~ x ) 



Pqg ~ Iqg 
Pgq = ^gq 



where 



l m = C F \ < —-l-x 
l q9 = (2a; -1) 

Igq = C F (2 - x) 

Igg = 2N c (-— + --2 + x{l-x)). 

V 1 — x x J 



(15) 



(16) 



These kernels need a suitable regularization to be well defined. Below we will analize the 
implicit regularization underlying eq. ()15|) . One observation is however almost immediate: 
the component P gg is not of the form given by eq. In general, therefore, in the singlet 
case, the generalization of eq. is given by 

P(x) = A(aO - ^(l - x) f 1 P 2 {z)dz (17) 

Jo 

and a probabilistic interpretation is more complex compared to the non-singlet case and 
has been discussed in the original literature [I]. 



4 Convolutions and Master Form of the Singlet 

Distributions are folded with the kernels and the result rearranged in order to simplify 
the structure of the equations. Since in the previous literature this is done in a rather 
involuted way p] we provide here a simplificaton, from which the equivalence of the various 
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is tne smipie relation 



1 d y ts i \ f 1 d vyf(y)- x f( x ) n \*r \ (iq\ 

f( x /v) = / + _ x )f\ x ) ( 18 ) 



2/(1-2/)+ ^ 2/ y-z 

in which, on the right hand side, regularity of both the first and the second term is explicit. 
For instance, the evolution equations become 

dq nn f dy yq(y) - xq(x) , nr ,, n w v f l dy 3 

— — — ^- = 2C F / + 2C F log(l - x) - / — (1 + z) q(y) + -C F q{x) 

dlog(Q 2 ) J y y — x Jx y 2 

+n f / 1 ^(, 2 + (l-,) 2 ),( 2/ ) 
Jx y x ' 

dy ,i 1 + (1 - Z) \ {y) + 2N C f - Xf ^ 9 {y) 



d\og(Q 2 ) Jx y z A y y _ x 

+2N C log(l - + 2N C f 1 ^ f- - 2 + z(l - *)) 4(2/) + ^(x) 

Jx y \z ) 2 



(19) 



with z = x/y. The same simplified form is obtained from the probabilistic version, having 
defined a suitable regularization of the edge point singularities in the integrals over the 
components 7/// in eq. (fTK|) . The canonical expressions of the kernels ([2*T]) . expressed in 
terms of "+" distributions, can also be rearranged to look like their equivalent probabilistic 
form by isolating the edge-point contributions hidden in their "+" distributions. We get 
the expressions 



and 

the other expressions remaining invariant. In appendix A we provide some technical 
details on the equivalence between the convolutions obtained using these kernels with the 
standard ones. 

A master form of the singlet (unpolarized) equation is obtained by a straightforward 
change of variable in the decreasing terms. We obtain 

t- = f —iw( x /y)(i(y)- I —i qq (y/x)q(x) 

clt Jx y Jo y 

t~ = [ —iggi x ly)- n f I % g (y/ x )g( x ) 

or Jx y Jo 

1 l-X—A rl dy 

~o lgg{y/ x )9(x) + —lgg(x/y)q(y) (22) 
z ja Jx y 
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01 tins aspect is leit in appendix id. ine ^reguiateuj transition prooauinties are tnen given 

by 

Wqg{x\y) = j gq (x/y) d( " V > X ^ V < 1 — — 
w (y\x) = 7 (v/x) %<^-A)%>0) 

6(y > x)6(y < 1 - A) 



w gg( x \y) = lggi. x /y)- 



Wqq(y\x) = (n nqg (y/x) - -Igg(y/X)j 
Wgq{y\x) = Jgq(x/y) 



y 

1 \ 0(y < x- A)% > 0) 



x 

6(y > x)6(y < 1 - A) 



y 

w gq( x \y) = °> 

(23) 

as one can easily deduct from the form of eq. (fT0|) . 

5 A Kramers- Moyal Expansion for the DGLAP Equa- 
tion 

Kramers-Moyal (KM) expansions of the master equations (backward or forward) are some- 
times useful in order to gain insight into the master equation itself, since they may provide 
a complementary view of the underlying dynamics. 

The expansion allows to get rid of the integral which characterizes the master equation, 
at the cost of introducing a differential operator of arbitrary order. For the approxima- 
tion to be useful, one has to stop the expansion after the first few orders. In many cases 
this turns out to be possibile. Examples of processes of this type are special Langevin 
processes and processes described by a Fokker-Planck operator. In these cases the prob- 
abilistic interpretation allows us to write down a fictituous lagrangean, a corresponding 
path integral and solve for the propagators using the Feynman-Kac fomula. For definitess 
we take the integral to cover all the real axis in the variable x' 

d f°° iii i i 

q(x, t) = / dx (w(x\x )q(x , r) — w(x \x)q(x 7 r)) dx . (24) 



dr 



— oo 



As we will see below, in the DGLAP case some modifications to the usual form of the 
KM expansion will appear. At this point we perform a KM expansion of the equation in 
the usual way. We make the substitutions in the master equation y — > x — y in the first 
term and y — ► x + y in the second term 

d f 00 

—q(x,T)= dy(w(x\x-y)q(x-y,T) -w(x + y\x)q(x,T)) , (25) 

<JT J — oo 

identically equal to 

— q(x,T) = j dy(w(x + y - y \x - y')q(x -y',r) - w(x + y \x)q(x,r)) , (26) 
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can oe reiateu using a tayior ^or rvivi j expansion 01 tne nrst term 



d r°° 00 (—y) n d n 

— q{x,r) = J^dy^Z n] dxn {w{x + y\x)q{x,r)) (27) 

where the n = term has canceled between the first and the second contribution coming 
from |26l) . The result can be written in the form 

Q 00 (—y) n d n 

frfa r ) = E\^^ Mx)q(x, r)) (28) 
where _ 

/oo 
dy(y-x) n w(y\x). (29) 
-00 

In the DGLAP case we need to amend the former derivation, due to the presence of 
boundaries (0 < x < 1) in the Bjorken variable x. For simplicity we will focus on the 
non-singlet case. We rewrite the master equation using the same change of variables used 
above 



d r 1 

q(x,r) = / dyw{x\y)q{y,r) - I dyw{y\x)q{x,r) 



Or 



where we have introduced the simplest form of the Moyal product 



a(x) 

dyw(x + y\x) * q(x,r) + / dyw(x + y\x)q(x, r), (30) 



2 



w(x + y\x) * q(x) = w(x + y\x)e 
and a(x) = x — 1. The expansion is of the form 



-y 



(^■ +? -)^,r) (31) 



' - / dyw(x + y\x)q(x,r) - V / dy — d x n (w(x + y\x)q{x, r)) (32) 

Ja(x) ^ Jo n! 



which can be reduced to a differential equation of arbitrary order using simple manipula- 
tions. We recall that the Fokker-Planck approximation is obtained stopping the expansion 
at the second order 

d 1 

Q-q(x, r) = a (x) - d x (a 1 (x)q(x)) + -d 2 x (a 2 (x)q(x, r)) (33) 

with 

a n (x) = / dyy n w(x + y,x) (34) 



being moments of the transition probability function w. Given the boundary conditions 
on the Bjorken variable x, even in the Fokker-Planck approximation, the Fokker-Planck 
version of the DGLAP equation is slightly more involved than Eq. (j3~3l) and the coefficients 
a n (x) need to be redefined. 

2 A note for noncommutative geometers: this simplified form is obtained for a dissipative dynamics 
when the p's of phase space are replaced by constants. Here we have only one variable: x 
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The probabilistic interpretation of the DGLAP equation motivates us to investigate the 
role of the Fokker-Planck (FP) approximation to the equation and its possible practical 
use. We should start by saying a word of caution regarding this expansion. 

In the context of a random walk, an all-order derivative expansion of the master 
equation can be arrested to the first few terms either if the conditions of Pawula's theorem 
are satisfied -in which case the FP approximation turns out to be exact- or if the transition 
probabilities show an exponential decay above a certain distance allowed to the random 
walk. Since the DGLAP kernels show only an algebraic decay in x, and there isn't any 
explicit scale in the kernel themselves, the expansion is questionable. However, from a 
formal viewpoint, it is still allowed. With these caveats in mind we proceed to investigate 
the features of this expansion. 

We redefine 



do(x) = / dyw(x + y\x)q(x, r) 

J a(x) 
ra(x) 

a n (x) = / dyy n w(x + y\x)q(x, r) 



o 



ra(x) 

d n (x) = / dyy n d x n (w(x + y\x)q(x, r)) n = l,2, ... (35) 



o 



For the first two terms (n = 1, 2) one can easily work out the relations 

di(x) = d x ai(x) — a(x)d x a(x)w(x + a(x)\x)q(x, r) 

a>2{x) = dla 2 (x) — 2a(x)(d x a(x)) 2 w(x + a(x)\x)q(x, r) — a(x) 2 d x a(x)d x (w(x + a(x)\x)q(x, r)) 
-a 2 (x)d x a(x)d x (w(x + y\x)q(x, r)) \ y=a ( x ) (36) 



Let's see what happens when we arrest the expansion ()32j) to the first 3 terms. The 
Fokker-Planck version of the equation is obtained by including in the approximation only 
d n with n = 0, 1, 2. 

The Fokker-Planck limit of the (non-singlet) equation is then given by 

d .... 1 

— q(x, t) = a (x) + a^x) - -a 2 (x) (37) 

which we rewrite explicitely as 

8 . . „ /85 3 13 10 12 nl fl-x\\ . . 

^4 + 8^-fa-T + ^r)^- r >- (38) 

A similar approach can be followed also for other cases, for which a probabilistic 
picture (a derivation of Collins-Qiu type) has not been established yet, such as for h\. We 
describe briefly how to proceed in this case. 



9 



lorm. ±ms is possiuie since trie suutraction terms can ue written as integrals 01 a positive 
function. A possibility is to choose the transition probabilities 

Wl [x\y] = -f-\-2]%> a ;)%<l) 

wma = v(r^-0 %> - x)%<o) 



(39) 



which reproduce the evolution equation for hi in master form 



dh\ 

— = dywi(x\y)hi{y,r) - dyw 2 {y\x)h l (x,T). (40) 
or Jo Jo 

The Kramers-Moyal expansion is derived as before, with some slight modifications. The 

result is obtained introducing an intermediate cutoff which is removed at the end. In this 

case we get 

— C ( U 2 I 3 6 I 21ogf — - X) h (x r) 

dr V 3 3a; 3 x 2 x \ x J J 

Notice that compared to the standard Fokker-Planck approximation, the boundary 
now generates a term on the left-hand-side of the equation proportional to q(x) which is 
absent in eq. (|33|) . This and higher order approximations to the DGLAP equation can be 
studied systematically both analytically and numerically and it is possible to assess the 
validity of the approximation [15j. 



7 Helicities to LO 

As we have mentioned above, an interesting version of the usual DGLAP equation involves 
the helicity distributions. 

We start introducing [2] the DGLAP kernels for fixed helicites P ++ (z) = (P(z) + 
AP(z))/2 and P+-(z) = (P(z) - AP(z))/2 which will be used below. P(z) denotes 
(generically) the unpolarized kernels, while the AP(z) are the longitudinally polarized 
ones. These definitions, throughout the paper, are meant to be expanded up to NLO, the 
order at which our numerical analysis holds. 

The equations, in the helicity basis, are 



2tt y y 

+Pf + (-) ® g + (y) + Pf_(~) ® 9-(v)l 

y y 
= ^(P+-(-) ® ff+(y) + p++(-) ® q-(v) 
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Figure 2: LO kernels (qq and qg) in the helicity basis. 



,x s 



,x s 



+P q + 9 -{-)®g + {y) + P q A{-)®g-{y)l 

dg+{x) _ as 
2n' 



dt 



,x s 



y 



+Pi 9 + (-)®9 + (y) + n 9 -(-)®g4y)), 
y y 



dg-jx) 
dt 



a 



,x s 



,x s 



^{Pn-)® q+ {y) + Pl\{-)®q-{y) 

+P 9 + 9 -{-)®g + {y) + P 9 + 9 + {-)®9-{y))- 
y y 



(41) 



The non-singlet (valence) analogue of this equation is also easy to write down 

dj ^ 1 = ^( p +-( X -) ® *+Av) + ^) ® (42) 

where the = q± — q± are the valence components of fixed helicites. The kernels in 
this basis are given by 



p(0) 


. p(0) 


_ p(0) 


p(0) 


p(0) 
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p(0) 


= n/x 2 
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= n f (x 
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(1-x). 



1 — x — x 2 

x 



PoS(l-x) 
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Figure 3: NLO kernels (qq and qg) in the helicity basis. 

P 5 ( 9 V = N c (3x+^-3-x 2 ^j (43) 

Taking linear combinations of these equations (adding and subtracting) , one recovers 
the usual evolutions for unpolarized q(x) and longitudinally polarized Aq(x) distributions. 
We recall that the unpolarized distributions, the polarized and the transversely polarized 
qr{x) are related by 



q(x) = q+ (x) + g_(x) = q +T (x) + q^ T {x) 
Aq{x) = g + (x) — q~(x) (44) 

at any Q of the evolution and, in particular, at the boundary of the evolution. 

Similar definition have been introduced for the gluon sector with G±(x) denoting 
the fixed helicities of the gluon distributions with Ag(x) = g+{x) — g~(x) and g(x) = 
g + (x)+g_(x) being the corresponding longitudinal asymmetry and the unpolarized density 
respectively. 

8 Summary of Positivity Arguments 

Let's recapitulate here the basic arguments |2| that are brought forward in order to prove 
the positivity of the evolution to NLO. 
If 

\AP(z)\ < P{z),z < 1 (45) 

then both kernels P ++ (z) and P + -(z) are positive as far as z < 1. 

The singular contributions at z — 1, which appear as subtraction terms in the evolution 
and which could, in principle, alter positivity, appear only in diagonal form, which means 
that they are only contained in -P++, multiplied by the single functions q+(x) or q~{x) 
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Figure 4: Nonsinglet kernels in the helicity basis. 
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Figure 5: LO kernels (gq and gg) in the helicity basis. 
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Figure 6: NLO kernels (gq and gg) in the helicity basis. 



dq- (x) 



dt 

dg+(x) 

dt 



as pqg i x * 



®g+(y) + ... 



(46) 
(47) 
(48) 



Let's focus just on the equation for q + (jlfij) . Rewriting the diagonal contribution as a 
master equation 



dq + (x) 
dt 



dx' (w ++ (x\x')q + (x' , r) — w ++ (x'\x)q + (x, r)) dx' + .. 



in terms of a transition probability 



a. 



w++(x\y) = —P ++ (x/y)6(y > x) 
In 



(49) 



(50) 



which can be easily established to be positive, as we are going to show rigorously below, as 
far as all the remanining terms (the ellipses) are positive. We have performed a detailed 
numerical analysis to show the positivity of the contributions at x — 1. 

This last condition is also clearly satisfied, since the 5(1 — z) contributions appear only 
in P ++ and are diagonal in the helicity of the various flavours (q, g). For a rigorous proof 
of the positivity of the solutions of master equations we proceed as follows. 

Let q(x, t) be a positive distribution for r < r c and let us assume that it vanishes at 
t = r c , after which it turns negative. We also assume that the evolution of q(x, r) is of 
the form (0) with positive transition probabilities w(x\y) and w(y\x). Notice that since 
the function is continuous together with its first derivative and decreasing, continuity of 
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requires tnat at r = t c 



dq(x) 
dt 



I 



dx'w(x\x')q(x', t) 



(51) 



which is positive, and we have a contradiction. We can picture the evolution in r of 
these functions as a family of curves getting support to smaller and smaller x- values as r 
grows and being almost vanishing at intermediate and large x values. We should mention 
that this proof does not require a complete probabilistic picture of the evolution, but just 
the positivity of the bulk part of the kernels, the positivity of the edge point subtractions 
and their diagonality in flavour. From Figs. (J2J and JSJ) it is also evident that the leading 
order kernels are positive, together with the qg and qq (Fig. |j3J)) sectors. 

The edge point contributions, generating the "subtraction terms" in the master equa- 
tions for the "++" components of the kernels are positive, as is illustrated in 3 Tables 
included in Appendix A. There we have organized these terms in the form ~ CS(1 — x) 
with 



with A and B being numerical coefficients depending on the number of flavours included 
in the kernels. Notice that the subtraction terms are always of the form (|52|) . with the (di- 
verging) logarithmic contribution (~ Jq dz/(l — z)) regulated by a cutoff. This divergence 
in the convolution cancels when these terms are combined with the divergence at x = 1 
of the first term of the master equation for all the relevant components containing "+" 
distributions. It is crucial, however, to establish positivity of the evolution of the helicities 
that the boundary conditions on the evolution \Aq(x,Ql)\ < q(x,Ql) be satisfied. Initial 
conditions have this special property, in most of models, and the proof of positivity of all 
the distributions therefore holds at any Q. 

As we move to NLO, the pattern gets more complicated. In fact, from a numerical 
check, one can see that some NLO kernels turn to be negative, including the unpolarized 
kernels and the helicity kernels, while others (Fig. (J3J)) are positive. One can also notice 
the presence of a crossing of several helicity components in the gq and gg sectors (see Fig. 
(|7|) at larger x values, while in the small-x region some components turn negative (Fig. 
Q). There is no compelling proof of positivity, in this case, either than that coming from 
a direct numerical analysis. 

9 NLO Numerical Tests using Recursion Relations 

We have seen that master forms of evolution equations, for evolutions of all kinds, when 
found, can be used to establish positivity of the evolution itself. 

The requirements have been spelled out above and can be summarized in the following 
points: 1) diagonality of the decreasing terms, 2) initial positivity of the distributions, 3) 
positivity of the remaining (non diagonal) kernels. As we have also seen, some of these 
conditions are not satisfied by the NLO evolution. 

In order then to proceed with a numerical test of the inequality we have decided to 
work directly in x-space, using a specific ansatz which summarizes the NLO evolution in 
a rather compact form. 

This ansatz, which we will illustrate below, reduces the evolution equations to a suit- 
able set of recursion relations [HI, El- 



C 



log(l - A)A + B 



(52) 



15 




Figure 7: The crossing of NLO kernels (gq and gg) in the helicity basis. 

In this ansatz, the NLO expansion of the distributions in the DGLAP equation is 
generically given by 



f< r*\ ^ A n (x) ( a(Q 2 )\ 2 ~ B n (x) n ( «(Q 2 ) 



n=0 



nl 



a(Q, 



n=0 



TV. 



(53) 



where we assume a short hand matrix notation for all the convolution products. Notice 
that f(x,Q 2 ) stands for a vector having as components all the helicities of the various 
flavours (q±,G±). The ansatz implies a tower of recursion relations once the running 
coupling is kept into account (THl EH] 



A n+l {x) = --P^ ® A n {x) 
Po 



(54) 



to leading order and 

B n+ i{x) = -BJx) 



4/?o 



^A n+1 (x)) - -^-P {1) ® A n (x) - |-P(°) ® B n (x) 
j ^Po Po 



~P {1) ® A n (x) - ® B n (x) 

Airpo Po 



(55) 



to NLO, relations which are solved with the initial condition B (x) = 0. The initial 
conditions for the coefficients A (x) and B (x) are specified with q(x, Qq) as a leading 
order ansatz for the initial distribution 



Aq(x) = 5(1 - x) ® q(x, Ql) = g (x) 



(56) 
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A (x) + a B (x) = q (x). 



(57) 



Again, other boundary choices are possible for A (x) and B (x) as far as (|57]l is 
fullfilled. 

Once the expansion is established, any linear combination of distributions has still to 
satisfy the same expansion, for appropriate combinations of the scale-invariant coefficients 
A n (x) given above. 

The numerical implementation of this algorithm, especially to NLO, is not so straight- 
forward and requires some effort. Among its advantages, however, is its remarkable speed. 



10 Initial conditions and Results 

As input distributions in the unpolarized case, we have used the parametrized formulas 
of Ref. [TT], that are calculated to NLO in the MS scheme at a scale Ql = 0.40 GeV 2 

x(u-u)(x,Ql) = 0.632x a43 (l-x) 3 - 09 (l + 18.2s) 

x(d-d)(x,Q 2 ) = 0.624(1 - xf- Q x(u -u)(x,Ql) 

x(d-u)(x,Ql) = 0.20x a43 (l -x) 12 - 4 (l - 13.3v^ + 60.0x) 

x(u + d)(x,Ql) = 1.24x a20 (l -x) 8 - 5 (l -2.3v^ + 5.7x) 

xg(x,Q 2 ) = 20.80x L6 (l-x) 4 - 1 (58) 

and xqi(x, Ql) = xqj(x, Ql) = for qi = s, c, b, t. 

Following [12], we have related the unpolarized input distributions to the longitudinally 
polarized ones by 



xAu(x,Ql) = 1.019x a52 (l -x)°- 12 xu{x,Q 2 ) 

xAd(x,Q 2 ) = -0.669x 0A3 xd(x,Ql) 

xAu(x,Ql) = -0.272x°- 38 xu(x,Ql) 

xAd(x, Ql) = xAu(x, Ql) 

xAg(x,Ql) = lA19x 1A3 (l-x)°- 15 xg(x,Ql) (59) 

and xAqi(x, Ql) = xAqj(x, Ql) = for qi = s, c, b, t. 

We show in Fig. |H1 results for the evolution of the u + distribution at the initial scale 
(0.632 GeV) and at two final scales, 100 GeV and 200 GeV respectively. The peaks at the 
various scales get lowered and become more pronounced toward the smaller x region as 
Q increases. In Fig. show an apparent steeper growth at small-x compared to u + . 

For the d distributions the situation is inverted, with d~ growing steeper compared to d + 
(Figs. [TT] and [12] respectively). This apparent behaviour is resolved in Figs. (fTOf and (fT3|) 
from which it is evident that both plus and minus components converge, at very small-x 
values, toward the same limit. 

The components s,c,t and b (Figs. rnH2*Tj) have been generated radiatively from van- 
ishing initial conditions for final evolution scales of 100 and 200 GeV (s, c, b) and 200 GeV 
(t). 
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larger x values. ±ne uistnuution tor tne top wuaiK yi) lias oeen inciuueu tor completeness. 
Given the smaller evolution interval the helicity distributions for heavier generations are 
suppressed compared to those of lighter flavours. Gluon helicities (Figs. [22J [231) are a l so 
enhanced at small-x, and show a similar growth. The fine difference between the quark u 
and d distributions are shown in Figs. fTUl and [TBI 

Finally in Figs. EUESlandEElwe plot simultaneously longitudinally polarized, unpolar- 
ized and helicity distributions for up quarks, down quarks and gluons at an intermediate 
factorization scale of 100 GeV, relevant for experiments at RHIC. Notice that while u + 
and u~ are positive and their difference (Au) is also positive, for down quarks the two 
helicity components are positive while their difference (Ad) is negative. Gluons, in the 
model studied here, have a positive longitudinal polarization, and their helicity compo- 
nents are also positive. The positive and negative gluon helicities are plotted in two 
separate figures, Figs. [22]and[2Hl while their difference, Ag(x) is shown in Fig. [23 One 
can observe, at least in this model, a crossing at small-x in this distribution. 

We conclude that, at least for this set of boundary conditions, positivity of all the 
components holds to NLO, as expected. 

11 Conclusions 

We have discussed in detail some of the main features of the probabilistic approach to 
the DGLAP evolution in the helicity basis. Numerical results for the evolution of all the 
helicities have been provided, using a special algorithm, based in x-space. We have also 
illustrated some of the essential differences between the standard distributional form of 
the kernels and their probabilistic version, claryfying some issues connected to their reg- 
ularization. Then we have turned to the probabilistic picture, stressing on the connection 
between the random walk approach to parton diffusion in x-space and the master form 
of the DGLAP equation. The link between the two descriptions has been discussed es- 
pecially in the context of the Kramers-Moyal expansion. A Fokker-Planck approximation 
to the expansion has also been presented which may turn useful for the study of formal 
properties of the probabilistic evolution. We have also seen that positivity of the helicity 
distributions, to NLO, requires a numerical analysis, as already hinted in [2j. Our study 
also validates the use of a very fast evolution algorithm, alternative to other standard 
algorithms based on Mellin algorithms, whose advantage is especially in the analysis of 
the evolution of nonforward parton distributions, as we will show elsewhere. 
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A Edge Point Positivity 

We report below 3 tables illustrating the (positive) numerical values of the contributions 
coming from the subtraction terms in the NLO kernels. Coefficients A and B refer to the 
subtraction terms — log(l — A) A + B as explained in the section above. 
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iz.i My 


4 


10.9569 


10.6924 


5 


9.3836 


9.2109 


6 


7.8103 


7.7249 


N f 


A 


B 


3 


12.5302 


12.1739 


4 


10.9569 


10.6924 


5 


9.3836 


9.2109 


6 


7.8103 


7.7249 


Nf 


A 


B 


3 


48.4555 


27.3912 


4 


45.7889 


24.0579 


5 


43.1222 


20.7245 


6 


40.4555 


17.3912 



Table 1. Coefficients A and B for P, 



(i) 



NS,++ 



Table 2. Coefficients A and B as in Table 1 for P. 



(i) 



qq,++ 



Table 3. Coefficients A and B as in Table 1 for P, 



(i) 



99,++ 



B Regularizations 

The "+" plus form of the kernels and all the other forms introduced before, obtained by 
separating the contributions from the edge-point (x — 1) from those coming from the 
bulk (0 < x < 1) are all equivalent, as we are going to show, with the understanding 
that a linear (unique) cutoff is used to regulate the divergences both at x=0 and at x=l. 
We focus here on the two possible sources of singularity, i.e. on P qq and on the P gg 
contributions, which require some attention. Let's start from the P qq case. We recall that 
"+" plus distributions are defined as 



1 



6(1 - x - A) 



— - d(l - x) / 

(l-a;) + l-x v ' Jo 1-z 

with A being a cutoff for the edge-point contribution. 



(60) 



We will be using the relations 

r = i 

Jx V 



1 dz 



o 1 — z 



I ~r~ log(l-x) 

Jx 1 — Z 

I —f(y)g(x/y) = [ —f(x/y)g(y). 
Jx y Jx y 



V Jx y 

Using the expressions above it is easy to obtain 



(61) 



1 



®f(x) 



dy 



-q(y) 



(1 - x)+ Jx y (1 - x/y). 

Jx U 1 — X / U Jx U JO Z 



y i - x/y Jx y 

1 dyyf(y) - xf(x) 



x y y-x 



+ \og(l-x)f(x) 



(62) 
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regularization, we neeu special care. 111 me prouaoinstic version 01 trie Kernel, rue nanunng 
of P qq (g) q is rather straightforward 



P qq ®q{x) = C F ['^(-^-l-x/y^qiy) 



(63) 



and using eqs. (ffiTJ) we easily obtain 

P qq ®q(x) = 2C F \-2C F \og(l - x)q(x) 

Jx y y — x 

-C F I' ^ (1 + s/y) g(y) + ~C F q{x). (64) 
Now consider the convolution P gg ® g(x) in the Collins-Qiu form. We get 

p 99 » 9 (x) = 2Ca /;| t -l_ 9W 

+2CU /' * (*/y(l - x/y) + ^- - 2) 9 (») - sto /* £ + ^ ~ ) 
Jx y \ x/y J Jo C a \z \ — zJ 

-~g(x) £ dz 2C A (z(l -z)-2)- n f g(x) £ d Z l - (f + (1 - zf) . (65) 

There are some terms in the expression above that require some care. The appropriate 
regularization is 

r 1 dz r 1 dz T ... r 1 dz r x ~ A dz 
/ - + / 1 > ' A = / — + / 1 • 66 

Jo Z Jo 1 — Z J A Z Jo 1 — Z 

Observe also that 

"i dz /-i-A dz 



, , , -log A. (67) 

I A Z Jo 1 — Z 

Notice that in this regularization the singularity of 1/z at z = is traded for a singularity 
at z=l in 1/(1 — z). It is then rather straightforward to show that 

Pgg ® #(x) = 2C A / h 2CU log(l - x)g(x) 

J x y y-x 

x/y(l - x/y) + -±- - 2^j + ^g(x). (68) 



+ 2cj ld -y( 

Jx y \ 



A final comment is due for the form of the sum rule 



dz{z - i) 7ra = (69) 



that we need to check with the regularization given above. 
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anu 1/^1 — z) integration terms anu use 

" l dz r 1 dz 



o z Jo 1 — z 



1 dz /-i-A 



a : Jo 1 — 2 
log A + log A = 



to eliminate the singularities at the boundaries x = 0, 1 and verify eq. (jEHj) . 
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Figure 8: Evolution of u + versus x at various Q values. 
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Figure 10: Small- a; behaviour of ^ at 100 GeV 
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Figure 11: Evolution of d + versus x at various Q values. 
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Figure 12: Evolution of d versus x at various Q values. 
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Figure 13: Small-x behaviour of d ± at 100 GeV. 
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Figure 14: Evolution of s + versus x at various Q values. 
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Figure 15: Evolution of s~ versus x at various Q values. 
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Figure 16: Evolution of c + versus x at various Q values. 
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Figure 17: Evolution of c~ versus x at various Q values. 



0.12 



0.1 h 



0.08 



0.06 



0.04 



0.02 



0=0.632 GeV 
Q=100 GeV 
Q=200 GeV 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x 

Figure 18: Evolution of b + versus x at various Q values. 
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Figure 19: Evolution of b versus x at various Q values. 



C Kernels in the helicity basis 

The expression of the kernels in the helicity basis given below are obtained combining the 
NLO computations of [T7 | ITfil PTS] 
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Figure 20: Evolution of t + versus x at various Q values. 
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Figure 21: Evolution oft - versus x at various Q values. 
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Figure 22: Evolution of g + versus x at various Q values, 
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Figure 23: Evolution of g~ versus x at various Q values. 
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Figure 24: Various kinds of distributions of quark up at Q = 100 GeV. 

+Nl (2(70 - 3x)x - 3tt 2 (1 - x + 3x 2 - x 3 )) 
+2N c T f (x(37 + x(23x - 57)) - 23)]} 

+ S.2C F T f (l - 3x) -N%U- 13x + j | logx 



+ 



2x 



1 — x + 3x — x ) > log x 



{2N 2 1 
(x 3 - 3x 2 + x - lH logxlog(l - x) 

\ N 2 C {\ + x(2 + 2x + x 3 )) 



x(l + x) 



So(x) 



(82) 



References 

[1] J.C. Collins and J. Qiu, Phys.Rev.D39 1398,1989. 

[2] C. Bourrely, E. Leader, O.V. Teryaev, |hep-ph/9803238| 

[3] C. Bourrely, J. Soffer, O.V. Teryaev, Phys.Lett.B420 375,1998. 

[4] D. Boyanovsky, H.J. De Vega, D.S. Lee, S.Y. Wang, H.L. Yu, Phys.Rev.D65 
045014,2002. 



[5] B. E. Baaquie, C. Coriano and S. Marakani, cond-mat/0211489 



[6] E. Eichten (Fermilab), I. Hinchliffe, Kenneth D. Lane and C. Quigg, 
Rev.Mod.Phys.56, 579, 1984. 



[7] M.Gluck, E.Reya and A.Vogt, Eur.Phys.J.C 5 (1998) 461 



32 



0.4 - 




xd(x) 
x Ad(x) 
x (f (x) 
xd (x) 



0.3 - 



0.2 \ 



0.1 - 



- 



-0.1 - 



-0.2 







0.1 



0.2 



0.3 



0.4 



0.5 



0.6 



0.7 



0.8 



0.9 



x 



Figure 25: Various kinds of distributions of quark down at Q = 100 GeV. 

[8] G. Rossi, Phys.Rev.D29:852, 1984. 

[9] J.H. Da Luz Vieira, J.K. Storrow, Z.Phys.C51 241,1991. 
[10] L.E. Gordon and G. P. Ramsey Phys.Rev.D58 094017, 1998. 
[11] M.Gluck, E.Reya and A.Vogt, Eur.Phys.J.C 5, 461, 1998. 

[12] M.Gluck, E.Reya, M.Stratmann and W.Vogelsang, Phys.Rev.D 63 094005, 2001. 

[13] A. Cafarella, C. Coriano and M. Guzzi, hep-ph/0209149, Proceedings of the Intl. 
Workshop on Nonlinear Physics, ed. B. Prinari et. al., World Scientific 2003 (in 
press). 

[14] O.Martin, A.Schafer, M.Stratmann and W.Vogelsang, Phys.Rev.D 57 1117502, 1998. 
[15] A. Cafarella, C. Coriano, in preparation. 

[16] G. Curci, W. Furmanski and R. Petronzio, Nucl.Phys.B175:27, 1980 
[17] R. Mertig and W. L. Van Neerven, Z.Phys.C70 637, 1996. 
[18] W. Vogelsang, Nucl.Phys.B475, 47, 1996. 



33 



Q=100 GeV 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x 

Figure 26: Various kinds of gluon distributions at Q = 100 GeV. 
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Figure 27: Small-x behaviour of Ag at 100 GeV. 
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